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Spatial dependencies of the pair potential and the local density of states near the surfaces 
of 6^2 _ y 2 -wave superconductors are studied theoretically. The calculation is based on the t-J 
model within a mean-field theory with Gutzwiller approximation. Various types of surface 
geometries are considered. Similar to our result in the extended Hubbard model, it is found 
that the formation of zero-energy states strongly depends on the surface geometry. In addition 
to this feature, the zero-energy states give peak splitting for the (110) surfaces when the super- 
exchange interaction J is large. This is due to the induced s-wave component near the surface. 
The present result explains the microscopic origin of the spontaneous time-reversal symmetry 
breaking at the surfaces of high-T c superconductors. 
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The t-J modeft is believed to explain the phase diagram of high-T c materials including the so- 
called pseudogap above T c and thus is a realistic model for high-T c super conductor sJ§ Although 
there are many works regarding this model, the quasiparticle properties in nonuniform systems 
including surfaces or interfaces are not clear at this stage except for those in the case with a 
vortexll) Since d x 2_ y 2-wave pair potential is the most promising symmetry in the t-J model for 
actual doping concentrations, we expect an interference effect of the quasiparticle due to the sign- 
change of the pair potential^ at the surface. 

Recent theoretical works within quasiclassical approximation revealed a novel interference effect 
peculiar to the unconventional superconductivity .i'i'Bi>i>0'0) One of the interesting interference 
effects is the appearance of zero-energy states (ZES) at the surface, which is due to the sign-change 
of pair potential when the quasiparticles are reflected at the surface.^ This localized ZES manifests 
itself as a zero-bias conductance peak (ZBCP) in tunneling experiments.!'© A recent experiment 
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using a well-oriented (110) surface detects a ZBCP.0) However this ZBCP is often not detected 
experimentally, which indicates that the presence of ZBCP strongly depends on the surface quality. 
One of the main purposes of the present letter is to clarify the origin of this dependence on the 
subtle difference of the surface in the t-J model. 

Another interesting feature is the splitting of ZBCP which was observed in recent experi- 
ments. This splitting is considered to be due to the coexistence of a s-wave pair potential 
induced near the surface. Till date, theories@lll§) have assumed an additional attractive potential 
which favors the s-wave component. For example, Fogelstrom et aZ.© suggested electron-phonon 
interaction as an origin of the s-wave component. However, in this letter, we will show that the 
t-J model automatically induces the s-wave component near the surface and thus gives a natural 
explanation of the splitting of ZBCP. We discuss the doping and J/t dependence of this splitting. 

Although the previous theories!'i'0!'!'00'll!@00'lil) clarified some important properties of 
high-T c superconductors, they ignored several distinctive features characteristic to high-T c mate- 
rials: i) short coherence length, which invalidates the quasiclassical approximation and ii) strong 
correlation. In order to deal with the effect of short coherence length, Tanuma et al& developed 
local density of states (LDOS) theory based on the extended Hubbard model beyond quasiclas- 
sical approximation. They clarified that LDOS of the d x 2 _ y 2-wave superconductor is sensitive to 
the atomic structures near the surface. However, the strong correlation effect is not sufficiently 
considered in their approach. 

To consider this problem, we develop a theory of the t-J model near the surface.© By using this 
model, we can naturally treat the d-wave superconductivity with short coherent length and strong 
correlation. The merit of studying the t-J model is that we can systematically investigate the 
doping or Fermi surface dependence which can be directly compared with the experimental results. 
In this letter, the LDOS of quasiparticles is obtained based on the self-consistently determined pair 
potential. Since it is difficult to treat the constraint in the t-J model analytically, we apply the 
Gutzwiller approximation© The validity of this method was verified by comparing the variational 
energies of the bulk states with those obtained in the variational Monte Carlo method© After 
the Gutzwiller approximation, the spatial dependence of the pair potential is determined self- 
consistently within the mean-field approximation as in our previous work regarding vortex coreJi^ 

The t-J model is given as& 

W = -* E (cUa + H.c.) + J ]T S i -S j 

i,cr 

where J and Si are the super-exchange interaction and the spin-1/2 operator at the i-th site, 
respectively. Here stands for the summation over the nearest-neighbor pairs. The operator 
c«r = Ci a PQ with the Gutzwiller projection operator Pq = n^l — ra^n^) excludes the double 



2 



occupancy. We use the Gutzwiller approximation!^ in which the effect of the projection is taken 
into account as statistical weights. The expectation values are estimated as 

( C L C J>) = 9t(4 a c ja ) , 
(S l -S j )=g s (S i -S j ) , (2) 

where (• • •) and (• • -)o represent the expectation values in terms of a Gutzwiller- type variational 
wave function Pg\&) and a BCS wave function, |<&), respectively. The renormalized coefficients 
gt = 25/ (1 + 5) and g s = 4/(1 + <5) 2 with the hole concentration 5(= 1 — n) are determined from the 
probabilities of involved configurations. EE EH Using this approximation, eq. (1) can be transformed 
into 

Ti-eS = -t c s ( c L c jv + H - c + J efi Si ■ Sj, (3) 



t c s = git, J e & = g s J- (4) 

In this letter, we consider various types of boundaries as shown in Fig. 1. The index m in Fig. 1(d) 
denotes the period of zigzag structures, and we will show that the LDOS strongly depends on this 
period. The case of m = [m = 1] corresponds to a flat (100) [(110)] surface shown in Fig. 1(a) 
[(b)]. In the following, we discuss the cases with m = 0, 1 and 2. 

For this model, we perform a mean-field approximation with site-dependent pair potential Ay 
and Hartree-Fock parameter 

Aij = -^Jcsic^Cjd, i ija = {c\ a c ja ). (5) 

Here we have assumed = £yj = £y. For simplicity £y and fi are fixed to the values £o an d 
jiQ determined in the bulk without boundaries. We use periodic boundary conditions in the y- 
direction and open boundary conditions in the x-direction. Furthermore, we assume that Ay is 
translationally invariant in the tangential direction along the surface. Thus, the unit cell is Nl x 1, 
with Nl being the number of sites in the x-direction. After Fourier transformation, the mean-field 
Hamiltonian becomes 

*j - ST ( n\ fu ^ n ( u \\ ( ^ {ky ^ ^ ky) \ ( Cttlfo) \ (a) 



with 
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Hij(k y ) = - £[(*cff + jJeB^o)(Sij±i + e^ lk y a 5i. j±m )} - noSij, (7) 
± 4 

&ij( k y) = Yl'Aij.Aj:\ + \j,ye Tikya 5 itj ± m ], (8) 
± 
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where Cj a (k y ) is the Fourier transformed form of Cj a with respect to the surface direction and 
j is now the site number in the x-direction (J = 1, • • • ,Nl). The wave vector k y changes from 
—ir/a to ir /a. In the above, a is the lattice constant and we use Nl = 300. The above Hamil- 
tonian is diagonalized by Bogoliubov transformationS© given as C^(k y ) = Y^i>lt(ky)Uj v , and 
Cii(-ky) = J2ult(^y)^N L +i v The s P a ti a l dependence of the pair potential is determined self- 
consistently as 

A j±md , y = \ J cff Uj^k+^mA 1 ~ f(Eu(ky))}e ±ik ^, (9) 

ky,V 

%fcl jVB = \ Jeff J2 Uj^k+jilA 1 ~ f@v(ky))}. (10) 

ky,V 

where Aj±m,j,y and Aj±i t j )X are the pair potentials along the y-axis and x-axis directions, respec- 
tively. The above f{E v (k y )) denotes the Fermi distribution function. 

We solve the effective Hamiltonian in Eq. (6) by numerical diagonalization and carry out an 
iteration until the pair potentials A%,j,x and Ajj i2/ are determined self-consistently. The obtained 
equations are decomposed into real and imaginary parts as 

ArjXj,) = Re(Ay i:I .( 2/ ))/Ao, 

\j,x(y) = Im (A-y, X (y))/A . (11) 

Figure 2(b) shows the calculated results of the spatial dependence of the pair potential for the flat 
(110) surface for J/t = 0.4 and 5 = 0.15. For this geometry, A^j x ^ y -j and Ajj x f y ) can be regarded as 
d-wave and extended s-wave components of the pair potential, respectively. The quantity A^j^y-j 
is suppressed at the surface and increases monotonically as we approach the middle of the lattice. 
This behavior is consistent with the quasiclassical theory© The extended s-wave component is 
induced near the surface, whose magnitude is about 30% relative to the bulk d-wave component. 
We find an atomic-scale spatial oscillation of the s-wave component, which is completely neglected 
in the quasiclassical approximation. 

Using the self-consistently determined pair potential, we calculate the LDOS at every site. In 
order to compare our theory with scanning tunneling microscopy (STM) experiments, we assume 
that the STM tip is metallic with a flat density of states, and that the tunneling probability is 
finite only for the nearest site from the tip. The LDOS at i-th site is given as, 

dwpi(o;)sech 2 (— — — ), (12) 
-co 2ksT 

p i (u) = --1m'£G?(k,u) 

= 2Y,Y.\ U iA 2 ^~Eu{k)) (13) 

k v 

where Gf'(k,u>) is the Fourier component of the retarded Green's function with energy uj. In the 
actual STM experiments, since the magnitude of the transparency between the tip and surface is 
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small, the tunneling conductance converges to the normalized LDOS 



m 




f , f s u2 ^ + 2A ^ 



(14) 



at low temperatures^ where p i: s{uj) denotes the LDOS in the superconducting state and Pn{^) 
denotes the LDOS in the normal state. In this letter, pn(uj) is obtained from the LDOS at the 
Nl /2-th. site away from the boundary. 

Figure 2(c) shows the calculated LDOS for various sites near the flat (110) surface. This zero- 
energy peak (ZEP) is the manifestation of ZES, which are formed due to the sign-change of the 
d x 2_ y 2-wave pair potential. A remarkable difference between the present results and those based on 
the quasiclassical theory is seen in the oscillatory behaviors of the LDOS© This oscillation can 
be regarded as the Friedel oscillation, the period of which is the inverse of the Fermi momentum. 
Furthermore, we find that the ZEP of the LDOS is split into two at all sites near the surface (see 
Fig. 2(c)). The splitting of the ZBCP in tunneling spectra is also obtained in the quasiclassical 
approximationBIll' Its origin is the s-wave component induced near the surface which blocks the 
motion of quasiparticles near the (110) surface.i'0) In effect, the splitting is not visible when J/t is 
decreased and the amplitude of the induced imaginary component is small. On the other hand, when 
5 is increased, the splitting becomes larger. This is because Arj^^) is reduced in magnitude with 
the increase of 6, while Aj j is insensitive to the change of 5. Recently, experimental observations 
of the peak splitting are reported in the tunneling spectroscopy of high-T c super conductor s.0'0) 
This indicates that the s-wave component is induced near the sample surface and the spontaneous 
time-reversal symmetry breaking actually takes place in high-T c superconductors. 

We also study the flat (100) surface shown in Fig. 1(a). In this case, Ai t j jX ( y \ is not induced near 
the surface. Since the quasiparticles do not feel the sign-change of the pair potential at the (100) 
surface, the ZEP do not appear. 

Next, we discuss the case of a 1 x 2 zigzag surface (see Fig. 3(a)). The obtained pair potential has 
complex spatial dependence as compared to the cases with (100) and (110) surfaces. The quantity 
A.i t j x f y \ is not induced as shown in Fig. 3(b). On the other hand, A R j x ^ has an oscillation and 
behaves like a d x 2_ y 2 pair potential far away from the boundary. The complex spatial dependencies 
of the pair potential reflect on the LDOS as an anomalous structure with many dips and peaks (see 
Fig. 3(c)). It should be noted that the ZES are not formed near the surface. A similar property 
was also obtained on the 1x2 zigzag surface in the extended Hubbard model.il) The reason for the 
absence of ZEP is as follows. The wave function of ZES spatially oscillates with the period of the 
inverse of the Fermi momentum as discussed in the (110) surface .III In the underdoped region, 
since the Fermi surface is nearly square, the period of the oscillation of the wave function is roughly 
coincident with 2a. Consequently, the node and the antinode appear alternatively. However, for 
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the 1x2 zigzag structure the phase of node and antinode does not coincide. This disappearance of 
the ZES can be regarded as an interference effect of the standing-wave, which cannot be explained 
by means of the quasiclassical theory. 

In this letter, we have investigated the LDOS near the surfaces of the d x 2_ y 2-wave superconductor 
based on the t-J model within the Gutzwiller approximation. In the present calculation, the 
non-local feature of the pair potential and the atomic-scale geometry of the surface are explicitly 
taken into account. The present result gives a microscopic basis for the spontaneous time-reversal 
symmetry breaking in d x 2_ y 2-wave superconductors where the s-wave component is induced as 
the d + is state near the surface. It is clarified that when the amplitude of the induced s-wave 
component is enhanced with the increase of the magnitude of J, the ZEP of LDOS is split into 
two. 
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Fig. 1. Schematic illustration corresponding tolxm zigzag surface : (a) a flat (100) surface (m = 0), (b) a fiat 
(110) surface (m = 1), (c) a 1 x 2 zigzag surface (m = 2) and (d) a 1 x m zigzag surface. 



Fig. 2. (a) Schematic illustration of a flat (110) surface, (b) the spatial dependence of the normalized pair potential 
and (c) the local density of states for J/t = 0.4 and S = 0.15 (Ao/t = 0.105). 



Fig. 3. (a) Schematic illustration of a 1 x 2 zigzag surface, (b) the spatial dependence of the normalized pair potential 
and (c) the local density of states for J/t — 0.2 and 5 = 0.1 (Ao/i = 0.054). 
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